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Abstract: Explicit equations are obtained to convert Cartesian coordinates to elliptic coordinates, based on which a function in 
elliptic coordinates can be readily mapped in physical space. Application to Kirchhoff vortex shows that its elliptical core 


induces two counter-rotating irrotational eddies. 
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1. Introduction 


elliptic coordinates (&, n) 


Figure 1. Coordinate lines for an elliptic coordinate system with a=1,b=0.5, 


elliptical coordinate system 


The (En) as a 


two-dimensional orthogonal coordinate system has many 
dynamical and engineering applications, such as Kirchhoff 
vortex [1], insect aerodynamics [2], hydrodynamic wave 


diffraction [3], and theoretical physics [4]. Its coordinate lines 
are confocal ellipses and hyperbolae and the transformation 
from elliptic to Cartesian coordinates is given by 


x = ccosh(é)cos(7) 
y =csinh(€)sin(7) 9 
Zg-b 


E20, 0<n<2z 


where a and b denote the semi-major and semi-minor axes of 
the ellipse and c is the elliptical eccentricity (Figure 1). 

In the meantime, explicit equations to transform from 
Cartesian to elliptic coordinates have not been found in the 
existing literature [5, 6, 7]. Such a conversion relation would 
be useful in mapping an elliptic-coordinate solution, for 
example the flow field of Kirchhoff vortex, in physical space. 


2. Cartesian to Elliptic Coordinates 


In order to invert the functional relation (1), we first 
eliminate ¢ and have 


2 2 


x V 
a ona C 
cos’(7) sin’ (77) 
which means curves of constant 77 are hyperbolae. The focus 
distance is c and the eccentricity is e= sec(77)- 


Let p = sin? (7), we have 
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which becomes 
cp’ +(x +y?-c’)p-y =0 (2) 
Then eliminating 77 from (1) we have 


2 2 


x + Y a 
cosh’ (E) sinh?” (é) 


It shows that curves of constant a are ellipses. The focus 


distance is c and the eccentricity is e = cosh” (£). 


Let g= —sinh? (£), we have 


2 


2 
l-q q 


which leads to 
Cg +(x? +y?-c’)q-y =0 (3) 


It is essentially the same as (2). Therefore (p,q) 


constitute the two roots of a quadratic equation. Since 
0<p<l, q<0,Wehave p> q, and the two roots are 


-B+ B° +4 y = 
Bae 2c’ a nS 


ae Ac? y? (4) 


in which B =x? +y =e 


From the definition of p we obtain 


h = aresin(,/ p ) (5) 


It has four cases depending on which quadrant the Cartesian 
point (x, y) is located, i.e., 


N=M)> x20,y 20 
N=T-h, x<0,y20 (6) 
N=2+M, x<0,y<0 
N=201-N, x>0,y<0 


Based on the definition of q, we can solve é from 


quadratic equation 
e“ +(4q—2)e** +1=0, 


which has two roots 
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e$ =1-2q+2Jq’ -q 


Since g <0, both roots are real and denoted as (¢,,€,). 


They clearly satisfy e’ eh =] , which leads to 


€,=-€ <0 . Because in elliptical coordinates only 


non-negative £ value is considered, we obtain 


§ =F In(1=29 +24" =4) (7) 


Eqs. (4-7) are explicit equations to derive elliptic 
coordinates from Cartesian grid. They can easily be realized 
via computation software such as Matlab. 


3. Application to Kirchhoff Vortex 


Kirchhoff vortex is a rotating elliptical region of uniform 
vorticity @ embedded in an irrotational ideal fluid [8]. It is 
the simplest example of non-smooth weak solutions to the 
Euler equations and has wide application in vortex dynamics 
[9-11]. It has a discontinuity of vorticity across its elliptical 
boundary 


2 y 
e 


S| x 


where a and b are semi-major and semi-minor axes, 
respectively. The boundary corresponds to contour 
d = arcsinh(b/c) , and the vortex rotates with constant 


angular velocity 


= ab a 
(a+by 


around the origin. From Act. 
streamfunction outside the core is 


159 of Lamb [1], the 


y= 72a +b) e* cos(29) + = wab ©) 


The instantaneous streamlines in a unsteady flow are given 
by the curves w= const. In a rotating frame with angular 


velocity © , Kirchhoff vortex looks steady and its 
streamfunction outside the elliptic core is related to the 
inertial-frame streamfunction (8) by 


m= y-5 En tén O 


From (9) it is straightforward to make a conformal mapping 
of steady streamfunction on a uniform mesh of elliptic 
coordinates (Figure 2). 
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Figure 2. Mapping of streamfunction (9) on a uniform mesh of (§,7), with a=1,b=0.5, @=10. 


In order to view the flow field in physical space, we use the 
conversion equations (4-7) to plot streamfunction (9) on a 
uniform Cartesian mesh in Figure 3, which shows the elliptical 
core of Kirchhoff vortex induces two irrotational eddies that 
rotate in opposite directions. 
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Figure 3. Mapping of streamfunction (9) on a uniform mesh of xy-plane, with 
a=1,b=0.5, @=10. The boundary of vortex core corresponds to § = 9.55, 


4. Conclusion 


This study obtains explicit equations that convert Cartesian 
coordinates to elliptic coordinates. The conversion relation 
can be easily realized with computer software and used to map 
a known function of elliptic-coordinates, such as the 
streamfunction of Kirchhoff vortex, on a uniform Cartesian 
mesh. 
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Matlab mfile scripts 


function [zeta,eta]=xy2el(x,y,a,b) 


% Conversion from Cartesian mesh (x,y) to elliptic coordinates (zeta,eta) 
% zeta (>=0), eta [0 2pi) 

% a -- semi-major axis 

% b -- semi-minor axis 

% c -- elliptic eccentricity 

% (x,y,zeta,eta) Same-size matrix 

% Author: Che Sun (CAS, 2016) 

% Source: http://eprint.las.ac.cn/abs/201611.00721 


c2=(a\2-bA2);c=2*c2; 

X2=x.A2; y2=y.A2; 

B=x2+y2-c2; 

del2=B.A2+2*c*y2; del=sqrt(del2); 
q=(-B+del)/c; q=sqrt(q); et0O=asin(q); 
eta=x; 


i=find(x>=0&y>=0); eta(i)=etO(i); 
i=find(x<0&y>=0); eta(i)=pi-etO(i); 
i=find(x<=0&y<0); eta(i)=pi+etO(i); 
i=find(x>0&y<0); eta(i)=2*pi-etO(i); 


p=(-B-del)/c; 
del2=p.\2-p; del=sqrt(del2); 
zeta=log(1-2*pt+2*del)/2;% only keep the positive root 


function [x,y]=el2xy(zeta,eta,a,b) 


% Conversion from elliptic coordiantes (zeta,eta) to Cartesian coordinates 
% zeta (>=0), eta [0 2pi) 

% a -- semi-major axis 

% b -- semi-minor axis 

% c -- elliptic eccentricity 

% (x,y,zeta,eta) Same-size matrix 

% Author: Che Sun (CAS, 2016) 

% Source: http://eprint.las.ac.cn/abs/201611.00721 


c2=a/2-b/2; c=sart(c2); 
x=c*cosh(zeta).*cos(eta); 
y=c*sinh(zeta).*sin(eta); 


